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Abstract 



I employ an exactly soluble toy model to investigate why unrenormal- 
ized perturbation theory works better than fully self-consistent approaches in 

describing the correlation-driven metal-insulator transition. 
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Theoretical works on the Mott transition [[!]] have mainly focused on the Hubbard model 
of strongly correlated electrons [Q] . Its main feature is the competition between itineracy and 
an on-site repulsion U. At half filling and zero temperature, the Hubbard model exhibits a 
metal-insulator transition at some critical value U c . At the moment, the only unified frame- 
work for describing the various features and phases associated with this Mott transition is 
provided by the dynamical mean-field theory ||, which becomes exact in the limit of infinite 
dimensions introduced by Metzner and Vollhardt [|J. However, the dynamical mean- field 
equations cannot be solved exactly, and various approximate numerical schemes have been 
employed to study their intricate structure. An important approximation for describing 
the half-filled Hubbard model and the Mott transition is the so-called iterated perturbation 
theory (IPT). This approximation succeeds to simultaneously capture the quasiparticle reso- 
nance and the incoherent lower and upper Hubbard bands. Partly based on the IPT, Zhang, 
Rozenberg, and Kotliar put forward a scenario for the Mott transition, in which the width of 
the quasiparticle resonance vanishes linearly in U — U c while its height remains unchanged, 
as the transition point U c is approached from below |5|]. This scenario was directly confirmed 
by the projective self-consistent technique of Moeller et al. ||. Yet, the IPT is based on 
second-order perturbation theory in U and might therefore be suspected to be inferior to 
perturbation-theoretical schemes that are either fully self-consistent or involve larger classes 
of diagrams. Self-consistent perturbation expansions have the advantage of conserving ex- 
act symmetries and their respective conservation laws ]/]]. Therefore, some authors have 
adopted a self-consistent approach from the outset |||§. Others, by contrast, have pointed 
out early on that fully self-consistent approaches seem to work worse than ordinary pertur- 
bation theory flOf . They relied on weak-coupling studies of the half-filled single-impurity 



Anderson model (TT|-|T3[. Meanwhile, there is a large body of numerical evidence, that, at 
least not too close to the Mott transition, IPT is quantitatively accurate P,|bl||. But it is not 



fully understood why leading-order perturbation theory works so well, whereas seemingly 
more elaborate self-consistent approaches notoriously fail to properly describe the insulating 
phase and their precursor effects. 



It is the purpose of this paper, to shed some light on this issue by considering an exactly 
solvable toy model that crudely captures some features of the Mott transition. This model 
is defined by the following Hamiltonian: 

H = U{n^ - + V£(f+c c + c+U (1) 

where h a = c+c CT . It describes a correlated c orbital hybridizing with an / orbital at zero 
energy. At half filling and zero temperature, particle-hole symmetry guarantees that the 
Green's function depends on three independent parameters only, 

e x = ~ [VU 2 + 64V 2 - VU 2 + 16V 2 ) , (3) 

e 2 = - (Vt/ 2 + 64V 2 + VU 2 + 16V 2 ) , (4) 

1 / f/ 2 - 32K 2 \ 

4 V v /(f/ 2 + 64V r2 )(t/ 2 + 16V r2 ); 
where 2 is a complex frequency and a-i — 1/2 — a\. At V — 0, the Green's function has two 
poles of equal weight at ±U/2 and thus represents an insulator. For V U, t\ ~ 6y 2 /f/, 
£2 — U/2+ 10V 2 /U, and ai ~ 18V 2 /?7 2 . The appearance of spectral weight at ±ei, close 
to the "Fermi level" at zero energy, is the best approximation to a metallic phase which 
is possible within our model. In the limit V — > 0, this spectral weight disappears, thus 
simulating a "metal-to- insulator transition." If we think of the two 5 function contributions 
at ±ei to the single-particle spectrum as representing the quasiparticle resonance (QPR), 
its total weight is given by Z = 2a x ~ 36V 2 /U 2 and its width by T* = 2e x ~ ZU/3. 

From Eqs. (@)-©, we obtain the noninteracting Green's function Go and the self-energy 
£(*) = Go-'W - G-'(z), 

- U2 ( 1 1 \ _ ^ 2 -2 

LU) " T \V^w + VTw) ~ Tz 2 -w 2 - {7} 

The first quantity is completely determined by V 2 , while the second one also depends on 
U 2 . The poles of the self-energy at ±3V are precursors of the 1/z pole of the insulating 



phase. Note, that close to the transition point, V <C U, these poles are located in the gap 
of the single-particle spectrum between the QPR and the "Hubbard bands" . The reason is 
that, while the width T* of the QPR collapses linearly in ZU, the self-energy poles move 
more slowly towards the Fermi level at u — 0, according to a y/ZU behavior. Note, that 
at least the dependences on Z are in agreement with the scenario of Ref. || on the Mott 
transition in the Hubbard model: There, T* ~ ZD, while the self-energy acquires poles at 
about ±\f~Z D (D is the half bandwidth of the bare band). 

We now consider the self-energy as a functional of the noninteracting Green's function. 
Eqs. (H) and (|7|) imply 

iE(f) = U^Goit)} 3 , (8) 

which means that the IPT approximation || becomes exact in our toy model. 

Next, we ask what we would obtain in leading-order self-consistent perturbation theory. 
The leading-order skeleton diagram for the self-energy is equivalent to Eq. (|8|), except that 
we have to replace the noninteracting Green's function by the fully renormalized one of Eq. 
(0). In the frequency domain, we obtain 

£(*) = U 2 (af 



1 1 

z - (2ei + e 2 ) z + (2e 1 + e 2 ) 

1 1 

z-(e 1 + 2e 2 ) z + {e x + 2e 2 ) 
1 1 

TUj 1 

U-3e 2 z + 3e 2 

As we approach the "insulating phase," V — * 0, only poles at ±3U/2 with residues U 2 /8 
survive. This approximation, therefore, fails to bring about the 1/z pole characteristic of the 
insulating phase. Its failure to account for the insulating phase becomes even more evident if 
we translate the result of Eq. (|9|) into a single-particle spectral function. For V = 0, we find a 
spurious 5 function contribution at the Fermi level, that carries 90% of the spectral weight. 
This problem is not specific to our toy model. In the context of the infinite-dimensional 
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Hubbard model, second-order self-consistent perturbation theory suffered from exactly the 
same flaws ||. There, too, it always predicts the existence of a QPR with a spectral weight 
that tends to be too large. 

Finally, we express the exact self-energy as a functional of the full Green's function 
(Q) with the goal of elucidating why its expansion in powers of U might break down as 
we approach the "transition point," V — > 0. The skeleton expansion of the self-energy 
is obtained by regarding G as a functional of Go and U, G[Gq,U], and by substituting 
the inverted functional Gq[G,U] into T,[Gq,U}. This requires the knowledge of G for an 
arbitrary Go and thus the solution of the most general impurity model, which is hopeless. 
Nevertheless, the essential element of the skeleton expansion is that the self-energy is viewed 
as a functional of G and U rather than Go and U. Since in our case, Go is parametrized by a 
single parameter, V 2 , straightforward perturbation theory consists in expanding £(V 2 , U) in 
powers of U. By contrast, renormalized perturbation theory corresponds to first trading V 2 , 
which completely determines Go, for one of the three parameters that determine the Green's 
function (0), say, ei, and afterwards expanding £(ei,£7) in powers of U. The procedure of 
going over from £(V 2 , U) to £(ei, U) corresponds to performing a Legendre transformation. 
Each of the Eqs. @-(§D can be solved uniquely for V 2 : 



V 2 = ^{l0e 1 + ^64e 2 1 + 9U 2 ), (10) 



V 2 = ^{l0e 2 v /GJrU9^). (11) 



, 9 - 40ai + 80af + 3 4ai - 1 a/9 - 8ax + 16a 2 

V 2 = - \ U 2 . (12) 

1024oi(l - 2a x ) v ' 

Depending on which of these equations we insert into Eq. (0), we obtain a different 
parametrization of S[G, U\. If renormalized perturbation theory were well-behaved all the 
way down to the transition point, V — > 0, all parametrizations of S[G, U] would also be 
well-behaved when expanded in powers of U. Yet, the expansion of S(ei,C/) in powers of 
U involves an expansion of the square root in Eq. ([To| ) in powers of U/ex, which converges 
only for U < 8ei/3. But since e\ — > as V — > 0, this expansion has zero radius of con- 
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vergence as one approaches the point V = 0, which is the analog of the Mott transition. 
More specifically, this expansion fails to converge for V 2 /U 2 < (5 + 4\/2)/G4 ~ 0.167, or 
for Z < 9(4 - >/2)/28 ~ 83%. This result indicates that the skeleton expansion of the 
self-energy is increasingly misbehaved as the Mott transition is approached, even though 
the actual range of convergence seems to be widely underestimated by our toy model. 

In summary, I have demonstrated within a simple toy model that neither second-order 
self-consistent perturbation theory nor a skeleton expansion of the self-energy to all orders in 
U can account for the Mott transition properly. While the former approach always predicts 
a metallic phase and overestimates the spectral weight of the quasiparticle resonance, the 
latter one fails to converge in the vicinity of the transition. By contrast, ordinary leading- 
order perturbation theory in U is well-behaved and, in our simple toy model, even turns out 
to be exact. 

I thank G. Kotliar for many stimulating discussions and the friendly working environ- 
ment. I would also like to thank the Deutsche Forschungsgemeinschaft for financial support. 
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